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2 New Three-Nucleon-Force Terms 

1 Introduction 

Recently it became possible to explain differences between data and predictions 
of modern nucleon-nucleon (NN) potentials for the total neutron-deuteron (nd) 
cross section and the minimum of the differential nd cross section at higher en- 
ergies by incorporating the 27r-exchange Tucson-Melbourne (TM) three-nucleon 
force (3NF) [y |y]. However, the puzzling discrepancy between theory and data 
for the low-energy analyzing power Q cannot be explained by that 3NF. Since 
it was shown in Q that this low-energy analyzing-power puzzle cannot be solved 
by reasonable changes in the NN potentials, one obviously needs new 3NF mech- 
anisms in order to resolve this long-standing mystery. 

The number of possible operators that can be used in constructing a 3NF is 
much larger than in the NN-force case, and it is not practicable to examine them 
all in order to see which are important for the low-energy analyzing powers and 
which are not. Rather, we need a systematic scheme that tells us which terms 
are likely candidates to be of importance and which are not. One such approach 
is chiral perturbation theory (xPT), which provides a power-counting scheme for 
the strength of the various 3NF terms. 

As we will explain below, xPT predicts in lowest non-vanishing order (beside 
the usual 27r-exchange terms) two terms of pion-range - short-range nature and 
three terms of short-range - short-range nature. Short range means, for example, 
that a meson heavier than a pion is exchanged between two of the three nucleons. 
Since the naive expectation is that pion-range - short-range 3NF terms are more 
important than the ones of short-range - short-range nature, in this paper we 
will deal only with the two pion-range - short-range 3NF terms, which we will 
include in our momentum-space calculations for the 3N continuum. Our interest 
will be concentrated on the low-energy vector analyzing powers. 

In Section ^ we review the history of the 27r-exchange 3NF and its effect in the 
3N continuum. In Section |3| we discuss those 3NF terms that include the exchange 
of mesons heavier than pions and that have been tested in the 3N continuum thus 
far. 

The two 3NF terms of pion-range - zero-range nature that are predicted in 
lowest non-vanishing order of xPT are introduced in Section 3. In this section 
we also explain how we adapt these terms, making them of finite range in order 
to be consistent with the traditional potentials that we use as our NN force, and 
also with the TM force that we use as our 27r-exchange 3NF in this paper. At the 
end of this section we discuss what conventional nuclear field-theory models say 
about these two 3NF terms. 

Our approach to solving the Faddeev equations for the 3N continuum is re- 
viewed briefly in Section g. 

The results for the nd elastic-scattering observables that incorporate the new 
3NF terms are presented and discussed in Section |6|. 

Finally we sum up and conclude in Section 0. 



In Appendix A. we present the partial-wave decomposition (PWD) for the 



two new 3NF terms. 
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Figure 1. The Fujita-Miyazawa 3NF. 



2 Short Review of 27r-Exchange Three-Nucleon Forces 

Three-nucleon forces have been part of nuclear physics for more than 40 years. 
Reahstic models began with the Fujita-Miyazawa (FM) force [^], which describes 
the exchange of two pions with an intermediate Z\-isobar as depicted in Fig. ||. 
Figure Q shows only that part of the 3NF for which nucleon 1 is the middle 
nucleon (i.e., the nucleon at which the virtual pion scattering takes place). We 
will call this configuration V4 a. The full 3NF is then given by 



.(1) 



Vi = v^'> + vr' + V, 



(2) 



^(3) 



(1) 



Scattering the pion in Fig. |l] from nucleon 1 via a virtual Z\-isobar is not 
the most general process that leads to a 27r-exchange 3NF. The many ways to 
accomplish this are indicated by the blob in Fig. ^. Several Ansatze have been used 
to derive the 27r-exchange 3NF up to now. A short overview is given in Table |2.1| . 
We will not comment here on the different techniques and the underlying physical 
ideas that lead to these 3NF models. For that we refer the interested reader to 



our recent paper |11]. Rather, we will concentrate here on the differences between 
these models in terms of operator form and the effects of these different models 
in 3N calculations. 

Let us begin with the types of operators contained in various 3NF models in 
use today. We will use the familiar language of the Tucson- Melbourne force |P], 
which labels four 27r-exchange operators by its parameters a, 6, c, and d. The 



^The use of the index "4" for a 3NF is common usage for three-nucleon calculations, where the 
indices "i"_"3" denote the three NN-pair forces. 
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Figure 2. The 27r-exchange 3NF. 



year 


3NF 


characteristic 


a' 


h 


c 


d 


1957 


Fujita-Miyazawa[ |] 


isobars 





-1.15 





-0.29 


1979 


Tucson- Melbourne 1^ 


current algebra 


-1.03 


-2.62 


1.03 


-0.60 


1983 


Brazil [0] 


chiral Lagrangian 
+ (current algebra) 


-1.05 


-2.29 


(1.05) 


-0.77 


1983 


Urbanai 


isobars with additional 
phenomenological 
medium-range term 





-1.20 





-0.30 


1993 


Texas |] 


chiral perturbation theory 


-1.87 


-3.82 





-1.12 


1996 


RuhrPotgO[ 


non-chiral Lagrangian 


-0.51 


-1.82 





-0.48 



Table 2.1. Various 27r-exchange 3NF models in use today. 
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operator form of these terms (neglecting all overall factors and form factors) is 
given by 
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where we have used the notation implicit in Figs. 11] and & 

The main difference between all these models (in terms of operator structure) 
is the c-term in the TM force. Let's have a closer look at the c-term and rewrite 
it (neglecting the isospin dependence in Eq. @)) as 
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Thus the c-term can be decomposed into a 27r-exchange term with the same 
operator structure as the a-term plus a short-range - 7r-range term (marked 
"S'-R" in Eq. (|6|) ) . The inclusion of the 27r-exchange part of the c-term leads to a 
redefinition of a as 



a — 2 m^ c 



(7) 



which essentially means a change of sign for a: a' ~ —a. Therefore the difference 
between the TM model and the other 27r-exchange 3NFs is the former's short- 



range - vr-range part of the c-term. Arguments developed in |11| using chiral 
symmetry show that this short-range - vr-range part of the c-term should be 
dropped. Doing this (and accordingly replacing a by a'), one gets a "corrected" 
TM force that we will call TM' in what follows. 

Since the EM force has included only the 2tt-A mechanism, it has only the b- 
and d-terms. Also, for the same reason, the values of b and d are roughly half the 
size of the corresponding TM-model values. 

The Brazil force [0] is very closely related to the TM model; however, it 
doesn't have the c-term. (That is, in it's original form, although a c-term was 
later included in order to agree with the TM model; this is indicated by the 
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brackets in Table |2.1| .) The values for the parameters o', b, and d are very close 
to the TM values. 

The Urbana model Q has only the b- and d-terms with values for b and d 
like those of the FM model, since it also incorporates only the Z\-mechanism. 
In addition it has a phenomenological medium-range term, whose inclusion was 
motivated by nuclear- matter calculations. It turns out that this term plays no 
important role in the 3N system. 

Like the Brazil force the Texas force B has the a-, b-, and d-terms, but with 
values for a', 6, and d obtained from a fit to tt — A^ scattering and substantially 
larger than those used in the Brazil or TM models. 



Finally, the RuhrPot model |1C] also employs the a'-, b-, and d-terms, but 
with significantly smaller values for the parameters than those used in the TM 
force. 

Consequently one does not expect major differences in the effects of these 
27r-exchange 3NFs in the 3N system and, indeed, no major differences have been 
found; (except for the Texas force) all of the above-mentioned 27r-exchange 3NF 
models have been tested in 3N bound-state and continuum calculations, either by 
the Bochum-Cracow group or the Pisa group or by both. Even the original TM 
force with the c-term gives results similar to all the other models at low energies, 
for reasons explained below. The only differences to be found are related to the 
different strengths of these models due to different values for the parameters, as 
well as due to different choices for the vrNN form factors. 

Let us now study the effects of the different terms individually. For the triton 
this has already been done in Ref. [^]. There it was found that the largest 
contribution to the binding energy (60%-70%) comes from the 6-term, whereas 
the a'-term can be neglected. The rest comes from the c- and d-terms, where 
the relative importance of these two terms depends strongly on the chosen NN 
interaction (Ref. |]12| used the RSC and AV14 potentials). 

In order to get the experimental value of the triton binding energy in con- 
junction with various NN force models, one can adjust the cut-off parameter A 
in the ttNN form factors of the TM and Brazil forces ||l3|| . The vrNN form factors 
of these models have the form 

Not only does the denominator in Eq. (8) suppress momenta large compared to 
the cut-off parameter A, but because F is normalized to one at Q^ = —m^, the 
numerator (and thus F at low momenta) also varies with A. This form factor 
therefore acts to a certain degree like a strength factor of the 3NF, which allows 
one to adjust the value of the triton binding energy. 

The different versions of the Urbana 3NF are fitted to give the correct value 
for the triton binding energy when used with one of the Argonne NN potentials. 

Most remarkable is the RuhrPot model, because the RuhrPot 27r-exchange 
3NF (together with the RuhrPot NN force) gives essentially the experimental 
value for the triton binding energy without having any adjustable parameters 
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Figure 3. Effects of the various terms of the TM 27r-exchange 3NF on the vector spin-correlation 
coefficient Cxx and the nucleon-to-deuteron tensor spin-transfer coefficient Ky ^ at Eiab = 3 
MeV. Predictions are: AV18 (solid line), AVf84-b (short-dashed line), AV18-fd (long-dashed 
line), AViS-fabd (dotted line) and AV18-|-abcd (long-short dashed line). The 3NFs are switched 
on only in the J^ = f /2 channels with jmax = 2. 



For the 3N continuum we have studied the effects of the individual 3NF terms 



{a',b,d), in conjunction with the AV18 np potential [15|, on elastic-scattering 
observables at Eiab = 3 MeV. The result is that the most important role in the 
27r-exchange 3NF is played by the 6-term. The other terms may have significant 
contributions to some (small) observables, but they are always smaller than the 
6-term contribution. As typical examples we depict the vector spin-correlation 
coefficient Cxx and the nucleon-to-deuteron tensor spin-transfer coefficient Ky ^ 
in Fig. |. 

This dominance of the 6-term explains why all 27r-exchange 3NFs, including 
the TM force with the c-term, give essentially the same results for the 3N con- 
tinuum after being fitted to the triton binding energy, even for those observables 
that do not scale with the triton binding energy. Of course, this picture might 
change if we go to higher energies. 

Interesting within this context is that many low-energy observables, especially 
in elastic scattering, show the just-mentioned scaling behavior with the triton 
binding energy [^ (i.e., all predictions for these observables using Hamiltonians 
with different NN potentials and a 3NF fitted to the triton binding energy agree 
with each other). We had a closer look at this scaling phenomenon and found (not 
unexpectedly) that the scaling observables are those that show a nonnegligible 
3NF effect only in J = 1/2"^ (J being the total 3-body angular momentum and 
n the parity), whereas non-scaling observables show 3NF effects also (or only) 
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Figure 4. Contributions of the TM 3NF in various channels for the scahng observable da; and 
the non-scaling observable Ay. The prediction without 3NF is the solid line (AV18 only). The 
TM 3NF is switched on for J^ = 1/2+ (short-dashed line), J^ = 1/2+ (long-dashed line), and 
J^ < 5/2 (dotted line), respectively. 



for other values of J^ . As typical examples we show Cxx and Ay in Fig. ^ 



3 Three-Nucleon Forces Including Heavier Mesons 



The TM model was extended to incorporate the exchange of p-mesons in |1(;], 
which leads to a vr-p and a p-p 3NF [17|. These forces were studied together with 



the TM vr-TT 3NF in the 3N elastic scattering and breakup process |18|. 

The result of this study was that the p-p 3NF has no visible effects, whereas 
the -n-p 3NF always has an effect opposite in sign to the effect of the vr-vr 3NF, 
but smaller. Moreover, it appears that the effect of the vr-p 3NF is more or less 
proportional to the effect of the vr-vr 3NF. 

That result is somewhat surprising. It means that replacing a vr with a p in 
the 3NF produces the same effects, though smaller and in the opposite direction. 
In other words, the effective physics in the vr-p 3NF is roughly the same as in the 
vr-vr 3NF. Given that the p and the vr interact very differently, this is a surprising 
conclusion. 

In order to understand this let's recall that the most important term in a vr-vr 
3NF is the 6-term. If we examine the operator structure of the vr-p 3NF (see, for 
example, the last paper of Ref. |l^), we find that this force also contains a term 
with the same operator structure as the 6-term of the vr-vr-exchange TM 3NF, 
but with different parameters and form factors and with the vr-mass in one of 
the meson propagators replaced by the p-mass. Moreover, we find that the vr-p 
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6-terni has the opposite sign to the vr-vr 6-term. So it seems plausible that the TM 
TT-p 3NF could be dominated by its 6-term, also; at least, that would explain the 
pattern found in the 3N scattering observables. (As a side remark, we note that 
the TT-p 3NF also has a d-term, which has the same sign as the vr-vr d-term.) 
Whether our supposition about the vr-/9 6-term is true or not, the finding in 



[18 1 about the effect of the vr-p TM 3NF on 3N scattering observables shows that 



the vr-/j 3NF of |17| cannot be expected to contribute to the explanation of any 
discrepancies between experiment and theory like the ^j^-puzzle. This vr-p 3NF 
just weakens the vr-vr 3NF, but does not lead to any new effects. In order to 
explain current puzzles we need new 3NF terms coming from other physics. 

4 The Texas Force and the Texas — Los Alamos Three-Nucleon Force 

In addition to the conventional 27r-exchange 3NF terms (the a'-, 6-, and d-terms 
in the TM language), xPT in first non- vanishing order predicts two terms of vr- 
range - zero-range nature (called di- and (i2-terms in the language of the Texas 
force) and three terms of zero-range - zero-range nature (called ei-, 62-, and 63- 
terms) P]. Since the zero-range - zero-range terms may be less important than 
the vr-range - zero-range terms (as the p-p terms in the TM force are much less 
important than the vr-p terms), we discuss here only the vr-range - zero-range 
terms, and defer until later the purely short-range terms. The form of the former 
terms in momentum space is given by 
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In a traditional Hamiltonian, of course, there are no zero-range forces. These 
are an artifact of xPT. A realistic force would contain short-range components 
from the exchange of heavy mesons. Indeed, one can construct a di-term from the 
exchange of w or o" mesons, whereas the (i2-term gets contributions from p and 
Al exchanges. We can envision a zero-range di-term as caused by the exchange 
of a very heavy isoscalar meson, and the corresponding d2 term by a very heavy 
isovector meson. Consequently we extend the zero-range forces to finite range 
by filling in meson propagators and adding form factors. Equation (0) is then 
modified to 

T/(i) ^1 3A ^ ^, ^ ^, F{Q' ) 2n - - 

di 9A .^ ,^,2. ^ ^ ^ ^ F{Q^) -. ^ .,,N 



10 
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Figure 5. vr-range - short-range 3NF terms. 
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Equation (O) is simply a heavy-meson propagator that is normalized to 1 at 
Q^ = 0; in Eq. (14) that propagator is multiplied by the product of two monopole 
form factors; and in Eq. (15) it is multiplied by the product of two dipole form 
factors. These terms are depicted in Fig. ^ 

The form factors in Eqs. ( [T^ ) and ( [Tsl ) are chosen in such a way that the 
resulting 3NF matrix elements at low momenta will not depend strongly on the 
value of the cut-off parameter Agr ■ In order to stay consistent with the TM' 27r- 
exchange 3NF, which we will use, we keep the same form for the ttNN form 
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factors used in the TM 3NF: 




nQ') = -TFTT^ (16) 



which is normahzed to 1 at the pion pole (Q^ = —rn^). Due to the —m^ factor 
in the numerator of Eq. (16), the size of a 3NF matrix element containing this 
form factor depends even for low momenta on the value of the cut-off parameter 
A, just as the TM 3NF does. 

Since the di- and d2-terms can be associated with the exchange of many 
different heavy mesons, we interpret them as effective forces subsuming the effects 
of all such mesons contributing to the respective term. Thus the exact value for 
the mass m^r in the propagator is not important; it just has to be roughly the 
right size. We choose rusr to be the w-meson mass. 

Another quantity that is unknown for the d\- and (i2-terms is their strength, 
since \^T cannot predict it. The only thing xPT can say about the strength of 
these terms is that they should be "natural". 

In order to see what naturalness means we rewrite the dimensionful coupling 
constants di and d2 in terms of dimensionless ones Q 

* - f-A <'«' 

where A = 1 GeV, fj^ = 92.4 MeV, and ci and C2 are dimensionless. For ci and 
C2 to be natural means that their value should be on the order of 1. In practical 
terms the absolute values of ci or C2 should typically (and roughly) be numbers 
comparable to 1,2,... (or the inverse); their signs, however, are not known. 

So in order to use the di- and d2-terms predicted by a xPT-based (i.e., the 
Texas) force, we have to make adaptations such as those described above. In order 
to differentiate between that zero-range force and the finite-range terms written 
down in Eqs. (^) and (12), we refer to the latter as the Texas - Los Alamos force. 



At this point we should mention that the di- and d2-teTuis are predicted not 
only by xPT, but by conventional nuclear field-theory models as well. For exam- 
ple, the di-term is contained in [^] as vr and a exchanges or vr and lu exchanges 
with an intermediate N(1440) resonance (Eqs. (3.2a) and (3.2b) in |]l^]), while the 
d2-term is included in pH] as a ir-p Kroll-Ruderman term (Eq. (2.13f) together 
with the last term, the "4", of Eq. (2.15a) in |l§). 

The d2-terni is also included in the meson-theoretical RuhrPot ir-p 3NF [1^] 
as Eq. (A7), among many other terms. Here the mechanism for the d2-teTin is a 
NNirp vertex at one of the three nucleons. 

However, these papers include many other 3NF terms, as well, that are of 
higher order in power counting (i.e., they are smaller) than the di and d2 terms. 



Reference |16] includes 126 terms, for example, and the (i2-term is part of one of 
them. After having written down these terms the authors conclude that the d2- 
part is the dominant one. This exhibits the advantage of a power-counting scheme 
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like the one that is part of xPT: one knows from the very beginning which terms 
should be important and which should not. With conventional nuclear field-theory 
models there is no easy way to know this beforehand. 

The di- and (i2-terms are therefore not new in the sense that they have never 
been written down before. But they are new in the sense that for the first time 
they have been identified as the leading-order terms of the pion-range - short- 
range 3NF and can be used (see below) in a realistic calculation of the 3N con- 
tinuum. 

5 Calculating the Three-Nucleon Continuum 



It became possible for the first time in |^0[ to include a realistic 3NF in a rig- 
orous calculation of the 3N continuum for energies above the deuteron breakup 
threshold. The algorithm used in pH] was later replaced by the more effective one 



developed in pi]. 



The Faddeev equation with a 3NF included reads 

T = tPcP + {1 + tGo)vi^\l + P)^ 



+ tPGoT + {l + tGom'{l + P)GoT (19) 



where T is the Faddeev amplitude for which Eq. (|19|) has to be solved, t is the 
2-body t-matrix, Go is the free 3-nucleon propagator, P is the sum of a cyclic 



^(1) 



and an anti-cyclic permutation (operator) for the three nucleons [22|, V4 has 
already been defined in Section y, and 4> is the incoming state composed of a free 
nucleon and a free deuteron. Once Eq. (|l9| ) is solved one gets the elastic transition 
amplitude via 

U = PGq^ + PT + Vi^\l + P)<l) + Vi^\l + P)GoT (20) 

and the transition amplitude for the breakup process via 

Uo = {l + P)T (21) 



The 3NF is built into Eq. (|l9[) in a perturbative way. Solving Eq. (|1^) by iteration 



not only gives the different orders in T, but the different orders in V^ as well. 
This will become important later on. Of course we always iterate Eq. (Il^) until 
we reach full convergence. 

In addition to using the new algorithm for the Faddeev equations, we have 
also replaced the PWD for the 3NF used in ||2^ by one developed later in |2^ ]. 
The reason is that the original PWD is numerically unstable for higher partial 
waves than those required in the early work. The PWD for the di- and d2-teTms 



is performed in Appendix A. For details on the numerics see Ref. y]. 

6 Results 

The first thing one realizes when dealing with the di- and (i2-terms is that these 
terms appear to be big. By big we mean that one cannot get convergence for 
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Figure 6. The elastic neutron analyzing power Ay at Eiab ~ 3 MeV. The solid line in both 
figures represents the prediction for the NN potential (AV18) alone. In the left-hand figure the 
long-dashed line is the prediction for AV18 plus the di-term 3NF with ci — —1. The short- 
dashed line is the prediction of AV18 plus the d2-terni 3NF with C2 = 0.5. (Note that the 
short-dashed line overlaps almost completely with the solid line, and is nearly invisible.) In the 
right-hand figure the long-dashed line is the prediction of AV18 plus the TM' 3NF with the 
cut-off parameter A = 5.215 m,r. For the short-dashed line on top of TM' the di- and d2-terms 
have been added with parameters ci = — 1 and C2 ~ 0.5, respectively. Form factors and other 
parameters have been chosen as described in Section W. The circles are nd data from Ref. |2q]. 



the iteration of Eq. (19) when ci and C2 are significantly larger than 1 and 0.5, 
respectively. Indeed, the convergence criterion (Eq. (5.15) of Ref . p^ - the differ- 
ence between the last and before-last result of the iteration of Eq. (|l9| ) , averaged 
in a specific way, has to become less than an e) cannot usually be fulfilled for 
e = 10"*^ (our usual value for e), but only for e = 10"'^ or larger. 

This slower (or lack of) convergence can easily be related to the size of these 
3NF terms. For example, if one slowly increases the value of ci from below 1 or 
C2 from below 0.5, one can examine how the convergence becomes worse around 
(or above) 1 or 0.5, respectively. 

One finds that the smallest e for which the convergence criterion can be ful- 
filled increases with increasing ci and C2 . If one increases the values for ci and C2 
further one finally reaches the point where convergence is totally lost. 

The above-described loss of convergence happens first for J = 1/2^, whose 
convergence is always slowest due to the presence of the bound state. However, 
even for larger values of ci and C2 where there is no longer convergence in the 
1/2+ channels, the 3NF contribution to the Faddeev amplitude T is still signifi- 
cantly smaller than the NN-force contribution. Thus naively one would expect no 
difficulty with the convergence. However, iterating Eq. (^) mixes the perturba- 
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tion series for the 3NF with the iteration series for the solution of T and therefore 
we lose convergence even for relatively small 3NFs. 

It might be that this behavior is an artifact of our way of parameterizing 
the new 3NF terms (e.g., of our choice for the form factors). Unfortunately, we 
didn't have the resources necessary to study this. Since we intend to study only 
the cardinal effects of these new 3NF terms, this problem is not particularly 
important. We simply confine ourselves to smaller values for ci and C2 for which 
we can present fully converged results. 

Even if the convergence problem turns out not to arise from our choice of 
parameterization, it is still not a problem of principle. We would just have to use 
a different algorithm for the incorporation of a 3NF into the Faddeev equations, 
either the one we used in the past ||2^ or, due to the advance of super-computers, 
just treat the 3NF in a straight-forward way (i.e., solve the Lippmann-Schwinger 
equation driven by V4 and then use the corresponding t-matrix in the Faddeev 
equations). 

All calculations presented in the following have been performed with a some- 
what reduced accuracy. That means that the maximally allowed angular mo- 
mentum for the two-body subsystem has been reduced to jmax = 2 everywhere 
except for the inner states during the calculation of the 3NF matrix elements (see 



Ref. 1 23 1 for details), which have been reduced to jmax = 3. With these restric- 
tions we are still within 2-3% of the results for a fully converged calculation at 3 
and 10 MeV and not much worse at 50 MeV. This is sufficient for our purpose. 

In the following we concentrate on the neutron analyzing power Ay. Of course, 
if we want to find a solution for the Ay-puzzle we also need to look at the deuteron 
analyzing power iTn. We don't do this here for several reasons. First, we do not 
intend to present a solution, but rather to see if a solution is possible. Second, 
there are no nd data for iTu, but only pd data. Since we cannot include the 
Coulomb force in our calculations it makes more sense to concentrate on Ay here. 
Third, the effects on iTn are always very similar to the effects on Ay] there are 
no significant differences. This suggests that if we are able to describe Ay we 
probably will also describe iTu. 

In Figs. ^^ we show the nucleon analyzing power Ay for elastic neutron- 
deuteron scattering. The negative sign for ci was chosen in order to get an 
enhancement in the maximum of Ay instead of a decrease. Although the en- 
hancement by di of the maximum of Ay at 3 MeV (depicted in Fig. y) is not 
terribly large and Ay calculated with the di-term is still far from the data, this is 
nevertheless a very promising result, since one can easily reach the experimental 
points with a larger value for ci (i.e., more negative than —1) that is still natural. 
We will return to this point. 

As can be seen from Fig. |6| the d2-tevm. has practically no effect on Ay at 3 
MeV. This is ideal since one can fix the value of ci at 3 MeV without worrying 
about effects from d2. 

One should note that the 27r-exchange 3NF (TM') already causes a visible in- 
crease in the maximum of Ay. However, plausible choices for the parameters in the 
27r-exchange 3NF will not be able to explain the Ay-puzzle, since these parame- 
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Figure 7. Same as Figure M but for Eiah ~ 10 MeV. The circles are nd data from Ref. 



ters are much more restricted than they are for the new 3NF terms. Our choice for 
the cut-off parameter of the ttNN form factors in the TM' 3NF (A = 5.215 itItj) 
probably makes our 27r-exchange 3NF somewhat too small. But even with the 
recommended choice of the Tucson-Melbourne group, one can only expect that 
the effect of the 27r-exchange 3NF would be larger by about 20%. This means 
that with a 27r-exchange 3NF one can come nowhere close to the experimental 
values. 

At 10 MeV the situation is similar to 3 MeV. One can see, however, that the 
d2-term has a small effect on Ay at this higher energy. 

At 50 MeV the situation has totally changed. Firstly, there is no Ay-puzzle at 
this energy (the predictions of all modern NN forces agree reasonably well with 
the data). However, it should be emphasized that the experimental situation is 
not as clear as at the lower energies, since the error bars are significantly larger. 
This means that there is still some room for small 3NF effects at this energy. 

Secondly, the form of Ay has changed. The maximum that we have seen 
at the two lower energies has become a minimum here, and there are two less- 
pronounced maxima, one at each side of the minimum. This is due to the fact that 
at this energy the mechanisms that build up the analyzing power have changed. 
We know [^ that at 50 MeV the P-waves have become less important and that 
D- waves play a role. 

If we now examine Fig. ^ we see that both the di- and the (i2-terms have 
large effects on Ay at this energy. This does not come as a surprise. We know 
that a 3NF always is of shorter range than the one-pion-exchange part of the NN 
force. Therefore one naively expects that 3NFs become more and more important 
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Figure 8. Same as Figure M but for Eiah = 50 MeV. The circles are nd data from Ref. 



with increasing energy. That this is indeed true for the 27r-exchange 3NF has been 
shown in Ref. Q and more recently in Refs. ||l[ and [^ for the total and differential 
elastic cross sections, respectively. Since the new 3NF terms d\ and ^2 are of even 
shorter range than the 27r-exchange 3NF, one expects that their effects grow even 
faster with increasing energy. Our calculations suggest that this is true. 

Of course, in this case we will face serious problems in our attempt to find a 
solution for the j4j^-puzzle via 3NFs of shorter range than the 27r-exchange 3NF. 
We will discuss this below. 

There is one other interesting aspect of Fig. |8|. If we look at the left maximum 
in the left-hand figure of Fig. |^ we see that the effects both of the d\- and d2-terms 
in that maximum are moderate. However, if we add these two terms on top of the 
27r-exchange 3NF this combined 3NF has very large effects on that maximum, 
as can be seen in the right-hand figure of Fig. ^, whereas the effect of the 27r- 
exchange 3NF alone is almost zero. This means that in this special case a very 
strong constructive interference develops between the various 3NF terms. It tells 
us that in dealing with 3NFs one has to be very careful if one excludes certain 
3NF terms, since even if their individual effects are small, their interference with 
other 3NF terms might lead to surprisingly large effects. 



Next we want to quantify our findings. For this purpose we list in Table 6.1 



the effects of the various 3NFs alone and together in the extrema of Ay. This 
table includes more calculations than we have shown in the figures. We included 
one calculation with c\ = —2 at 3 MeV. For this calculation we could only reach 
a convergence with e = 10~^. Therefore we did not repeat this calculation for the 
higher energies. 
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TM' - di 
+0.5d2 



3MeV 



10 MeV 



50 MeV 



50 MeV 





■Ay\max\Wo ) 


^y\m,ax\^'^'^-'^ ) 


^j;Uax (47.5°) 


^j/|mm(ll'^-5 j 


AV18* 
AV18 


0.04549 
0.04536 


0.1275 
0.1294 


0.1599 
0.1539 


-0.5066 
-0.5016 


+TM* 
+TM' 


0.04857(+6.8%) 
0.04825(+6.4%) 


0.1328(+4.2%) 
0.1420(+9.7%) 


0.1543(-3.6%) 
0.1472(-4.5%) 


-0.4953(-2.3%) 
-0.5066(+1.0% ) 


+di 
-di 

-2di 


0.04343(-4.4%) 

0.04778(+5.3%) 

0.05094(+12.3%) 


0.1231(-5.1%) 
0.1370(+5.9%) 


0.1589(+3.2%) 
0.1388(-10.9%) 


-0.4306(-16.5%) 
-0.5606(+11.8%) 


+0.5d2 
-0.5d2 


0.04521(-0.3%) 
0.04561(+0.6%) 


0.1277(-1.3%) 
0.1469(+13.5%) 


0.1578(+2.5%) 
0.1503(-2.4%) 


-0.5157(+2.8%) 
-0.4832(-3.8%) 



0.05106(+12.6%) 0.1497(+15.f 



0.04804 



-0.5367(+7.0%) 



Table 6.1. Theoretical predictions for Ay at the extrema for various energies. The cm. angles of 
the extrema are given in brackets in the second header line of the table. For calculations where 
a 3NF is added to AV18, the deviation of that result from the calculation with AV18 alone is 
given in per cent in brackets as well. All calculations were performed with jmax ~ 2 except the 
ones in the two lines marked with *, which were performed with jmax = 3. For simplicity we 
denoted different values for the dimensionless constants ci and C2 as multiplicative factors in 
front of di and d2 in the first column of the table (e.g., -2di means ci — —2). 



The first thing we see in Table |6.l| is that the effect of the elimination of the 
short-range part of the c-term in the original TM 3NF is negligible at 3 MeV, 
but already significant at 10 MeV. This means that older calculations with the 
original TM 3NF should be repeated using TM' for energies above the deuteron 
breakup threshold. 

Next we see that the effect of the di-term at 3 and 10 MeV is more or less 
linear with the value of ci, but not so at 50 MeV. On the other hand the effect of 
the (i2-term is not linear with the value of C2 at 3 MeV and is dramatically less 
so at 10 MeV, but becomes more or less linear at 50 MeV. 

A look at the last line of Table |6.l| reveals that at all energies there are some 
interference effects between the various 3NF terms, though they are strongest in 
the maximum at 50 MeV. Interestingly, we have a significant destructive inter- 
ference in the minimum of Ay at 50 MeV. 

Next we want to see if it is possible to find a combination of di and ^2 with 
which it is possible to come close to the experimental data for Ay. Of course, 
due to the nature of our calculations, any such combination that we might find 
can only give a rough estimate for the values of ci and C2, but that is all we 
want. If we can find such a combination we would have shown that a solution of 
the ^y-puzzle using these two 3NF terms of pion-range - short-range nature is 
possible. 

So let us start with 3 MeV. Here the situation is relatively simple, since we 
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can neglect the effect of the d2-term on Ay. Taking into account that we probably 
underestimate the effect of the 27r-exchange 3NF somewhat (as we stated above), 
we can extrapolate from Table |6 .11 that a value for ci of about —3 would probably 
be able to close the gap (~ 30%) between the data and the predictions of the NN 
potentials. 

If we move on to 10 MeV we see that (not taking into account d2 for the 
moment) a combination of TM'—3di (this notation means that we are using TM' 
and a di force with ci = —3) would probably overestimate the data a little bit, 
since the gap is also about 30% at this energy, but the 27r-exchange 3NF has a 
somewhat larger effect at 10 MeV than at 3 MeV, whereas the effect of di is more 
or less the same at both energies. However, at 10 MeV the (i2-term has a small but 
visible effect on Ay. So we can use the (i2-term to counterbalance the increased 
effect of the 27r-exchange 3NF. This leads us to a value for C2 that should be 
positive and small, perhaps 0.5 or 1. Thus a combination like TM' —3di + 0.5^2 
or TM'—Sdi + d2 would be able to bring the theoretical prediction for Ay close 
to the experimental data at 3 and 10 MeV. 

It might be, however, that the values we just gave for ci and C2 are even 
smaller in reality. The reason is that, as mentioned above, we used the ttNN 
form factor Eq. (16) for the di- and (i2-terms with the same small value for A as 
for the 27r-exchange part of the 3NF. This might cause us to underestimate the 
strength of di and ^2 somewhat. 



Now let's look at 50 MeV. We see immediately from Table |6j that the com- 
bination of 3NF terms mentioned above does not work at this energy. The effects 
would be far too large. If that is so, is there any way out of this dilemma? 

One might argue that 3NF terms of even shorter range than di and ^2 might 
become important at this higher energy. If that is the case these additional 3NF 
terms will provide additional parameters with which one might be able to describe 
the data at 50 MeV. As we mentioned in Section such terms exist: the so-called 
ei-, 62-, and e3-terms, which are of short-range - short-range nature. However, 
we believe that the inclusion of these terms will lead to serious problems for two 
(connected) reasons. 

We face the conceptual and practical problem that the short-range ei-, 62-, 
and es-terms are indistinguishable from very short-range (i.e., (^-function) parts 
of the 27r-exchange 3NF. Disentangling these terms is usually made more serious 
by the lack of consistency between the NN-force and the 3NF models we use. In 
addition, accumulating too many parameters (the values of ei-, 62-, and 63 are 
unknown) in a three-nucleon problem is self-defeating. 

This leads us immediately to a deeper, basic question: how high in energy 
do models for the nuclear force based on meson exchange make sense? Meson- 
exchange models have been extremely successful, even at much higher energies 
than we are considering here. Could it be that the meson-exchange picture - or 
rather the way we usually implement it - breaks down (i.e., becomes excessively 
complicated) at an energy as low as 50 MeV as soon as one starts to look into 
(admittedly small) details of the nuclear force? To answer this question one would 
have to investigate whether the short-range - short-range 3NF terms play any 
role in Ay (and other observables) at 50 MeV. Since this is such a basic question 
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it might be worthwhile to do so. 

However, there is one other possibihty. We have not yet mentioned the so- 



caUed 3NF Born terms |25, [Tlj. These terms are also predicted by xPT to be 
the same order in power counting as all the other 3NF terms mentioned here. 
These terms are model dependent in the sense that they depend on how the 
"off-shellness" of the NN potential has been defined, and they thus depend on 
the details of the NN potential that one uses [p8| , |ll[ . Usually these Born terms 
are neglected, on the one hand due to the complications they bring (such as 
nonlocality), and on the other hand because a subset of them are known to be 
small Q. One of the 16 Born terms is of spin-orbit nature [28|, however. This 
particular Born term might be important for Ay, since that observable is very 
sensitive to the spin-orbit force [^. The next step in exploring the ^y-puzzle 
should be to take into account (at least) this particular Born term. 

We have found that the effects of the di- and d2-tenQs on other observables in 
elastic nd scattering are usually smaller than their effects on Ay and iTn. Only 
for some tensor-polarization observables are the effects comparable to those on 
Ay and iTn. 

7 Summary and Conclusions 

In Section ^ we have reviewed the current state of affairs for the 27r-exchange 
3NF. Due to the elimination of the pion-range ~ short-range part of the c-term 
of the TM 3NF, the question of the operator form of the 27r-exchange 3NF is 
now settled. We showed that in low-energy elastic nd scattering the effect of the 
27r-exchange 3NF is dominated by the 6-term, as in the bound state. In addition 
we found that those observables that scale with the triton binding energy show 
3NF effects only for the channels for which the 3N system is bound (viz., those 
with J" = 1/2+). 

Next we commented in Section ^ on the results found so far using 3NFs that 
include the exchange of one or more mesons that are heavier than the pion. Due 
to the nature of these effects we suspect that the n-p TM 3NF is also dominated 
by its 6-term and therefore cannot explain the Ay-puzzle. 

In Section Q we introduced the Texas force, which is based on xPT. We 
explained which of the 3NF terms of the Texas force are of interest to us, and we 
extended the zero-range parts of these terms to finite range in order to connect 
with the traditional models we use for the NN force and the 27r-exchange 3NF. 
Thereafter we gave a meson-exchange interpretation to the new terms. Finally we 
commented on the appearance of those terms in calculations that are not based 
on xPT. 

A brief review of our approach to solving the Faddeev equations for the 3N 
continuum has been given in Section |5[ For the PWD of the new 3NF terms we 



refer to [Appendix A 



Finally we presented our results in Section ^ We studied the effects of the 
new 3NF terms on Ay at 3, 10, and 50 MeV. For the two lower energies we could 
find a combination of the di- and (i2-terms that (together with the NN force 
and the 27r-exchange 3NF) would be able to describe the Ay data. Although this 
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is only a qualitative finding, it is of considerable importance because we have 
developed for the first time a microscopic model of the nuclear force that has 
the potential to describe the low-energy vector-analyzing-power data. To make 
this model quantitative will involve considerable effort. We typically obtain the 
b and d parameters from fits to vr — A^ scattering, and in principle the di and d2 
parameters can be obtained from pion production in NN scattering, although the 
latter remains to be seen. If that proves impractical we would have to fit di and 
(^2 (and the cut-off parameter A in the ttNN form factor) to the triton binding 
energy, to Ay, and to other nd observables for which we have nd data. Given the 
complexity of three-nucleon calculations, this three-parameter search would be 
very time consuming. 

However, as discussed in Section ^, such a 3NF would not be able to describe 
the Ay data (and possibly data for other observables) at 50 MeV and higher. 
This might be due to the fact that other 3NF terms than those included in this 
paper become important at that energy. One such 3NF term could be one of 
the Born terms that is of spin-orbit type, and therefore is a candidate to be of 
importance for Ay. Up to now the Born terms have been largely neglected, but 
it seems obvious that at least the one just mentioned has to be included in 3N 
continuum calculations as a next step towards a solution of the Ay-puzzle. 

Other candidates for 3NF terms that might become important at 50 MeV are 
the ones of short-range - short-range type, which are predicted by xPT to be the 
same order as the di- and d2-terms and the 27r-exchange 3NF terms that have 
been included in this study. If it turns out that these terms do become important 
at 50 MeV we face serious problems with meson-exchange models. 

Unfortunately, another problem would also arise: we would have to fit many 
more terms, which on the one hand is very difficult for a 3N problem, and on 
the other hand we might loose any predictive power by accumulating too many 
parameters with increasing energy. Since these are important questions, it might 
be worthwhile to check first whether or not the short-range - short-range 3NF 
terms predicted by xPT in lowest non-vanishing order have any effect at 50 MeV. 

In this paper we believe that we have made an important step forward toward 
the solution of the long-standing Ay-puzzle by identifying new 3NF terms that 
have a significant effect on Ay even at low energies. We found these terms using 
a systematic approach for the classification of the 3NF terms (i.e., the power- 
counting scheme of xPT). In doing so new questions arose that require further 
testing of our models of nuclear forces. However, we have to leave the answers to 
these questions to future work. 
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For the partial- wave decomposition (PWD) of the new 3NF terms we will closely follow Ref. p3]. 
We will not repeat here the principles and ideas behind the PWD as developed in Ref. [eS but 
rather apply them to the new di and d2 terms. We will often refer to equations from Ref. P3|] in 
the form (xajpS]), where xx is the equation number in Ref. |23|. In addition, we refer to Ref. [P3| 
for the notation used here. 

A.l The di-Term 

Let us begin with the di-term. It is given in Eq. (nil) as 
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For the notation see Fig. H. As in Ref. |23] we will deal with the momentum - spin matrix 
elements and the isospin matrix elements separately (see Eq. (35[E3|)). 

We start with the momentum - spin matrix elements. In analogy to Eqs. (28|23|), (3ip3]) 
and (35123]) we split the momentum - spin dependent part of the di-term into two quasi- two- 
body operators. This leads to matrix elements 
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for the second term in Eq. (A.l). The momentum transfers Q and Q' are given by 

Q = P-Pi 
Q' = p' - P2 
For Mj' we have to decompose the operator 
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For the expansion of the angular dependence in the propagat or a nd form factor we use 
Eqs. (43-47pi) and (SgH). Putting this together with the result lK3) we get 
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Ne xt w e need to recouple the spherical harmonics that appear in the two matrix elements 
ofEq. ( IatI) : 
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Now we have to evaluate the two matrix elements in Eq. ( [A.10[ ): 
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Inserting these results into Eq. ( ^AW yields 
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Thus the final result for M^' , after evaluating the sums over m and mi in Eq. (A. 2), becomes 
fj,2 _ S{qi - q) 



M 



di 



Sjji 5mmi 5jji 5aai Sill 
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A useful test for the correctness of the result (A.14) is a comparison with the result for 



M;;'^ in PJ. The difference between M^^' 



and M^' is that the first matrix element has an 



operator ai ■ Q = in iai,Y^{Q)\ , whereas the latter matrix element has the operator 



Thus the difference between M^' and M^' , besides a differ- 



ent factor in front, is that in one case we have a rank-0 operator and in the other case the same 
operator, but with rank-1 this time. Of course the replacement of the rank-1 operator in the 
result for the d-term with the rank-0 operator in order to get the result for the di-term is not as 
straightforward as it might appear and has to be done with great care. Nevertheless, one finds 
that the result ( |A.1^ ) for A//;^ 

is consistent with the result (90p3|) for M^' . 

Besides M,' we will also need the matrix element 



^d'l =3 {P2q2a2j\ <Tl ■ Q' S2 ■ Q' JT 



F{Q'^) 



+ mi 



I ' ' ' \ 



(A.15) 



for the first term in Eq. (A.l). A-I^' can easily be obtained from M^' via the symmetry relation 



Mj^[ (pga,pigiai) =M^' {piqi 



ai,pqa) 



(A.16) 



The second momentum - spin matrix element occurring in the di-term is M^' , Eq. (A.3). 
For the moment we do not need to specify which choice of Eqs. (1l3-|lq) we want to make for 
the short-range operator 0{Q' ). That means we just use Eq. (46p3[) for the expansion of the 
angular dependence and do not yet specify how to determine H. With this we get 



^if = 2^ E ^r^K/pa) 3 {P2q2a2j\ If- [p'p2) \p q a'j) ^ 



(A.17) 



The matrix element in Eq. (A.16) can easily be evaluated as 



{P2q2a2j\ y ijiv'P'i) \p'q'a'j)^ = 5^,^^ Sp 



i-r 



and with this we get 
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(A.19) 



The symmetry relation to obtain M^ ' is given by 

^I'd^^iPQO'^Pili'^i) = M^f{piqiai,pqa) 

Let us now determine H for the different choices Eqs. (P_3l-P_5|) for 0{ 
according to Eq. (|l3|) as a propagator only we get instead of Eqs. (43-44| 



(A.20) 
Choosing 0(Q2) 
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and therefore H becomes for this case 
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The second case, Eq. (M), where the propagator is multipUed by monopole form factors, 
has aheady been calculated in Eqs. (43-47[p3|). Thus for this case one has 



' ppi 2{ppi) 



(A.23) 
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Finally we have the case where the propagator is multiplied by dipole form factors, Eq. (|l3) ; 
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This leads to 
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(A.25) 



The only pieces still missing for the PWD of the di-term are the isospin matrix elements. 
They are given by 



Id, = 2 {{tl/2)TMT\ ri ■ fa \{t'l/2)T' M't) ^ 

- S S , f 6^ x/fP / ^/2 1/2 f 1 r 1/2 1/2 t' 

id, = 2 {{tl/2)TMT\Ts ■ fi \{t'l/2)T'M^) 



(A.26) 



= 5' 



TT' OMr,,M' 



(-) 



t+t' + l 



6 Vtt 



1/2 1/2 i 1 j 1/2 1/2 t 
1/2 r i' J 1 1/2 1/2 1 



Id, for i = t' 
—Id, for f / f' 
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Now we can put all parts together in order to obtain the di matrix element. In order to do 
this we will use an obvious symbolic notation: 

(A.28) 



i{VrU, ) = Pi^2 {[M^[' B2«3 M^-;] Id, + [Mil' -P^^s Mif] Id,} P3«i 



A.2 The d2-Term 

The d2-term is given in Eq. (113) as 
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-^||.^x.-3.Q'a^2.Q'^^C?...(Q^)n.f2xf3 
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Ti-rsX fa (A.29) 



The isospin dependence is the same as in the d-term. The isospin matrix element for the 
da-term is therefore already given by Eq. (34|23]): Id2 ~ Id and Id2 = —Id- 

The d2-term does not fall into two parts as naturally as the di-term does. Nonetheless we ca n 
still split it into two quasi-two-body operators by rewriting (for the second term in Eq. (A.29)) 
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(A.30) 



in analogy to Eq. (77 P3|). Th erefore we have to calculate the following two matrix elements for 
the second term of Eq. (A. 29): 



fj,2 



F{Q^)\piqiai)2 



3{q-qi) 



^^Vf «{-'?/ >4' 

Sxxi 5m E C{jm IM - m, JM) C(jimi hlvh - mi,JiMi) 



5-i ■ Q 



X 2 



(pjml i-/2 y -^ Q |o"i. y ab^PP'^n 



1,-M 



(73 • ( 



+ mi 



F{Q^) buimi), (A.31) 



Mi'' = 3{P2q2a2\a^O{Q')\p'q'a')^ 

= ^^'^^ ,'; '^ ^ Syx, Sri, Yl C'O'a"^^ -^^Ma - ma, J2M2) C(j'm' /'A/' - m', J'A/') 



X 3 {P2J2m2\ a!^ 0{Q') |p'j'm') 



(A.32) 



The sum over fi from Eq. (A.3C) will be performed later on after we will have put together all 
the pieces. 

The first matrix element, M^' Eq. (A.31), occurred already in the d-term (with opposite 
sign) and is given by Eq. (90|23[): 
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(A.33) 



(Note that there are two misprints in Eq. (90|E3]): The phase in the second line of Eq. (90 1 
must read (— )'^^^^, and the phase (— )* in the third line has to be eliminated.). 
For the first term in Eq. (A.29) we need 

aixa3-g' = iV&J^-Q'\{ai,a:,}\Y\p'-p-2)Y 
= -i^]lfQ' f3,{ai,y^Q'}^}"° 
= iV2j^Q' J2^-r< {^.,Y\Q')} 
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which leads us to calculate the matrix element 

Mi^ = 3(p...a.hy2yf 0'{a.,y>V.)}'"'^^^F(Q-)|pVa')3 

(A.35) 

and Mf^ = -MJ'^. Eq. (91|||]) gives MJ'^ and Eqs. (92-93|2|]) the relation between Mf^ and 
M^' . Note that Eq. (93|E3|) has a misprint in the phase; the correct phase must be (— )"+"!. 
For the second matrix element, M^' in Eq. (A.32), we need 

Again, what w e have to insert for H depends on our choice for O. 



With Eq. (A.36) we get after performing the sums over m2 and m' in Eq. (A.32) 
Md2 = a a ^i2i' ^^2^ S12I' v6 27r Hi,(p2p ) 

/' J2 J'} C'(1M'^'^'.-^2^2) (A-37) 
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